Whole-exome sequencing reveals an association of rs112065068 in TGOLN2 gene with distant metastasis of non-small cell lung cancer.

Early prediction and prevention of recurring illness is critical for improving the survival rates of patients with non-small cell lung cancer (NSCLC). Previously, we demonstrated that the presence of premalignant epithelial changes in the small bronchi distant to the primary tumor is associated with NSCLC progression: isolated basal cell hyperplasia (iBCH) indicates a high risk of distant metastasis, BCH combined with squamous metaplasia (BCHSM) - a high risk of locoregional recurrence. Here, we aimed to identify germline single nucleotide variants (SNVs) and insertions and deletions (InDels) associated with distant metastasis and locoregional recurrence in cases with iBCH and BCHSM using whole-exome sequencing of 172 NSCLC patients. The rs112065068 of the TGOLN2 gene was identified only in iBCH patients and was associated with a high risk of distant metastasis (P < .001) and worse metastasis-free survival (HR = 4.19 (95 %CI 1.97-8.93); P < .001). This variant was validated in a group of 109 NSCLC patients using real-time PCR and Sanger sequencing analyses. To our knowledge, this study is the first to identify a germline variant associated with NSCLC distant metastasis.

CR1 variants contribute to FSGS susceptibility across multiple populations.

Focal segmental glomerulosclerosis (FSGS) is a common cause of nephrotic syndrome with an annual incidence in the United States in African-Americans compared to European-Americans of 24 cases and 5 cases per million, respectively. Among glomerular diseases in Europe and Latin-America, FSGS was the second most frequent diagnosis, and in Asia the fifth. We expand previous efforts in understanding genetics of FSGS by performing a case-control study involving ethnically-diverse groups FSGS cases (726) and a pool of controls (13,994), using panel sequencing of approximately 2,500 podocyte-expressed genes. Through rare variant association tests, we replicated known risk genes - KANK1, COL4A4, and APOL1. A novel significant association was observed for the gene encoding complement receptor 1 (CR1). High-risk rare variants in CR1 in the European-American cohort were commonly observed in Latin- and African-Americans. Therefore, a combined rare and common variant analysis was used to replicate the CR1 association in non-European populations. The CR1 risk variant, rs17047661, gives rise to the Sl1/Sl2 (R1601G) allele that was previously associated with protection against cerebral malaria. Pleiotropic effects of rs17047661 may explain the difference in allele frequencies across continental ancestries and suggest a possible role for genetically-driven alterations of adaptive immunity in the pathogenesis of FSGS.

Annotation of uORFs in the OMIM genes allows to reveal pathogenic variants in 5'UTRs.

An increasing number of studies emphasize the role of non-coding variants in the development of hereditary diseases. However, the interpretation of such variants in clinical genetic testing still remains a critical challenge due to poor knowledge of their pathogenicity mechanisms. It was previously shown that variants in 5'-untranslated regions (5'UTRs) can lead to hereditary diseases due to disruption of upstream open reading frames (uORFs). Here, we performed a manual annotation of upstream translation initiation sites (TISs) in human disease-associated genes from the OMIM database and revealed ∼4.7 thousand of TISs related to uORFs. We compared our TISs with the previous studies and provided a list of 'high confidence' uORFs. Using a luciferase assay, we experimentally validated the translation of uORFs in the ETFDH, PAX9, MAST1, HTT, TTN,GLI2 and COL2A1 genes, as well as existence of N-terminal CDS extension in the ZIC2 gene. Besides, we created a tool to annotate the effects of genetic variants located in uORFs. We revealed the variants from the HGMD and ClinVar databases that disrupt uORFs and thereby could lead to Mendelian disorders. We also showed that the distribution of uORFs-affecting variants differs between pathogenic and population variants. Finally, drawing on manually curated data, we developed a machine-learning algorithm that allows us to predict the TISs in other human genes.

Case report: Somatic mutations in microtubule dynamics-associated genes in patients with WNT-medulloblastoma tumors.

Medulloblastoma (MB) is the most common pediatric brain tumor which accounts for about 20% of all pediatric brain tumors and 63% of intracranial embryonal tumors. MB is considered to arise from precursor cell populations present during an early brain development. Most cases (~70%) of MB occur at the age of 1-4 and 5-9, but are also infrequently found in adults. Total annual frequency of pediatric tumors is about 5 cases per 1 million children. WNT-subtype of MB is characterized by a high probability of remission, with a long-term survival rate of about 90%. However, in some rare cases there may be increased metastatic activity, which dramatically reduces the likelihood of a favorable outcome. Here we report two cases of MB with a histological pattern consistent with desmoplastic/nodular (DP) and classic MB, and genetically classified as WNT-MB. Both cases showed putative causal somatic protein truncating mutations identified in microtubule-associated genes: ARID2, TUBB4A, and ANK3.

A 300-kb microduplication of 7q36.3 in a patient with triphalangeal thumb-polysyndactyly syndrome combined with congenital heart disease and optic disc coloboma: a case report.

BACKGROUND: Triphalangeal thumb-polysyndactyly syndrome (TPT-PS) is a rare well-defined autosomal dominant disorder characterized by long thumbs with three phalanges combined with pre- and postaxial polydactyly/syndactyly of limbs. By now, the syndrome has been reported in several large families from different ethnic backgrounds, with a high degree of inter- and intrafamilial variability. The genome locus responsible for TPT-PS has been mapped to the 7q36.3 region harboring a long-range sonic hedgehog (SHH) regulatory sequence (ZRS). Both single-nucleotide variants and complete duplications of ZRS were shown to cause TPT-PS and similar limb phenotypes. TPT-PS usually forms as isolated limb pathology not associated with additional malformations, in particular, with cardiovascular abnormalities. CASE PRESENTATION: Here we report on a rare Russian neonatal case of TPT-PS combined with severe congenital heart disease, namely double outlet right ventricle, and microphthalmia with optic disc coloboma. Pedigree analysis revealed TPT-PS of various expressivity in 10 family members throughout five generations, while the cardiac defect and the eye pathology were detected only in the proband. To extend the knowledge on genotype-phenotype spectrum of TPT-PS, the careful clinical and genomic analysis of the family was performed. High-resolution array-based comparative genomic hybridization (array-CGH) revealed a ~ 300 kb microduplication of 7q36.3 locus (arr[GRCh37] 7q36.3(156385810_156684811) × 3) that co-segregated with TPT-PS in the proband and her mother. The duplication encompassed three genes including LMBR1, the intron 5 of which is known to harbor ZRS. Based on whole-exome sequencing data, no additional pathogenic mutations or variants of uncertain clinical significance were found in morbid cardiac genes or genes associated with a microphthalmia/anophthalmia/coloboma spectrum of ocular malformations. CONCLUSIONS: The results support the previous data, indicating that complete ZRS duplication underlies TPT-PS, and suggest a broader phenotypic impact of the 7q36.3 microduplication. Potential involvement of the 7q36.3 microduplication in the patient's cardiac and eye malformations is discussed. However, the contribution of some additional genetic/epigenetic factors to the complex patient`s phenotype cannot be excluded entirely. Further comprehensive functional studies are needed to prove the possible involvement of the 7q36.3 locus in congenital heart disease and eye pathology.

Estimation of amyloid aggregate sizes with semi-denaturing detergent agarose gel electrophoresis and its limitations.

Semi-denaturing detergent agarose gel electrophoresis (SDD-AGE) was proposed by Vitaly V. Kushnirov in the Michael D. Ter-Avanesyan's laboratory as a method to compare sizes of amyloid aggregates. Currently, this method is widely used for amyloid investigation, but mostly as a qualitative approach. In this work, we assessed the possibilities and limitations of the quantitative analysis of amyloid aggregate size distribution using SDD-AGE results. For this purpose, we used aggregates of two well-characterized yeast amyloid-forming proteins, Sup35 and Rnq1, and developed a protocol to standardize image analysis and process the result. A detailed investigation of factors that may affect the results of SDD-AGE revealed that both the cell lysis method and electrophoresis conditions can substantially affect the estimation of aggregate size. Despite this, quantitative analysis of SDD-AGE results is possible when one needs to estimate and compare the size of aggregates on the same gel, or even in different experiments, if the experimental conditions are tightly controlled and additional standards are used.

No More Tears: Mining Sequencing Data for Novel Bt Cry Toxins with CryProcessor.

Bacillus thuringiensis (Bt) is a natural pathogen of insects and some other groups of invertebrates that produces three-domain Cry (3d-Cry) toxins, which are highly host-specific pesticidal proteins. These proteins represent the most commonly used bioinsecticides in the world and are used for commercial purposes on the market of insecticides, being convergent with the paradigm of sustainable growth and ecological development. Emerging resistance to known toxins in pests stresses the need to expand the list of known toxins to broaden the horizons of insecticidal approaches. For this purpose, we have elaborated a fast and user-friendly tool called CryProcessor, which allows productive and precise mining of 3d-Cry toxins. The only existing tool for mining Cry toxins, called a BtToxin_scanner, has significant limitations such as limited query size, lack of accuracy and an outdated database. In order to find a proper solution to these problems, we have developed a robust pipeline, capable of precise 3d-Cry toxin mining. The unique feature of the pipeline is the ability to search for Cry toxins sequences directly on assembly graphs, providing an opportunity to analyze raw sequencing data and overcoming the problem of fragmented assemblies. Moreover, CryProcessor is able to predict precisely the domain layout in arbitrary sequences, allowing the retrieval of sequences of definite domains beyond the bounds of a limited number of toxins presented in CryGetter. Our algorithm has shown efficiency in all its work modes and outperformed its analogues on large amounts of data. Here, we describe its main features and provide information on its benchmarking against existing analogues. CryProcessor is a novel, fast, convenient, open source (https://github.com/lab7arriam/cry_processor), platform-independent, and precise instrument with a console version and elaborated web interface (https://lab7.arriam.ru/tools/cry_processor). Its major merits could make it possible to carry out massive screening for novel 3d-Cry toxins and obtain sequences of specific domains for further comprehensive in silico experiments in constructing artificial toxins.

Census of halide-binding sites in protein structures.

MOTIVATION: Halides are negatively charged ions of halogens, forming fluorides (F-), chlorides (Cl-), bromides (Br-) and iodides (I-). These anions are quite reactive and interact both specifically and non-specifically with proteins. Despite their ubiquitous presence and important roles in protein function, little is known about the preferences of halides binding to proteins. To address this problem, we performed the analysis of halide-protein interactions, based on the entries in the Protein Data Bank. RESULTS: We have compiled a pipeline for the quick analysis of halide-binding sites in proteins using the available software. Our analysis revealed that all of halides are strongly attracted by the guanidinium moiety of arginine side chains, however, there are also certain preferences among halides for other partners. Furthermore, there is a certain preference for coordination numbers in the binding sites, with a correlation between coordination numbers and amino acid composition. This pipeline can be used as a tool for the analysis of specific halide-protein interactions and assist phasing experiments relying on halides as anomalous scatters. AVAILABILITY AND IMPLEMENTATION: All data described in this article can be reproduced via complied pipeline published at https://github.com/rostkick/Halide_sites/blob/master/README.md. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Phenome-wide functional dissection of pleiotropic effects highlights key molecular pathways for human complex traits.

Over the recent decades, genome-wide association studies (GWAS) have dramatically changed the understanding of human genetics. A recent genetic data release by UK Biobank (UKB) has allowed many researchers worldwide to have comprehensive look into the genetic architecture of thousands of human phenotypes. In this study, we used GWAS summary statistics derived from the UKB cohort to investigate functional mechanisms of pleiotropic effects across the human phenome. We find that highly pleiotropic variants often correspond to broadly expressed genes with ubiquitous functions, such as matrisome components and cell growth regulators; and tend to colocalize with tissue-shared eQTLs. At the same time, signaling pathway components are more prevalent among highly pleiotropic genes compared to regulatory proteins such as transcription factors. Our results suggest that protein-level pleiotropy mediated by ubiquitously expressed genes is the most prevalent mechanism of pleiotropic genetic effects across the human phenome.

Whole­exome sequencing in Russian children with non­type 1 diabetes mellitus reveals a wide spectrum of genetic variants in MODY­related and unrelated genes.

The present study reports on the frequency and the spectrum of genetic variants causative of monogenic diabetes in Russian children with non­type 1 diabetes mellitus. The present study included 60 unrelated Russian children with non­type 1 diabetes mellitus diagnosed before the age of 18 years. Genetic variants were screened using whole­exome sequencing (WES) in a panel of 35 genes causative of maturity onset diabetes of the young (MODY) and transient or permanent neonatal diabetes. Verification of the WES results was performed using PCR­direct sequencing. A total of 38 genetic variants were identified in 33 out of 60 patients (55%). The majority of patients (27/33, 81.8%) had variants in MODY­related genes: GCK (n=19), HNF1A (n=2), PAX4 (n=1), ABCC8 (n=1), KCNJ11 (n=1), GCK+HNF1A (n=1), GCK+BLK (n=1) and GCK+BLK+WFS1 (n=1). A total of 6 patients (6/33, 18.2%) had variants in MODY­unrelated genes: GATA6 (n=1), WFS1 (n=3), EIF2AK3 (n=1) and SLC19A2 (n=1). A total of 15 out of 38 variants were novel, including GCK, HNF1A, BLK, WFS1, EIF2AK3 and SLC19A2. To summarize, the present study demonstrates a high frequency and a wide spectrum of genetic variants causative of monogenic diabetes in Russian children with non­type 1 diabetes mellitus. The spectrum includes previously known and novel variants in MODY­related and unrelated genes, with multiple variants in a number of patients. The prevalence of GCK variants indicates that diagnostics of monogenic diabetes in Russian children may begin with testing for MODY2. However, the remaining variants are present at low frequencies in 9 different genes, altogether amounting to ~50% of the cases and highlighting the efficiency of using WES in non­GCK­MODY cases.

Whole-exome sequencing provides insights into monogenic disease prevalence in Northwest Russia.

BACKGROUND: Allele frequency data from large exome and genome aggregation projects such as the Genome Aggregation Database (gnomAD) are of ultimate importance to the interpretation of medical resequencing data. However, allele frequencies might significantly differ in poorly studied populations that are underrepresented in large-scale projects, such as the Russian population. METHODS: In this work, we leveraged our access to a large dataset of 694 exome samples to analyze genetic variation in the Northwest Russia. We compared the spectrum of genetic variants to the dbSNP build 151, and made estimates of ClinVar-based autosomal recessive (AR) disease allele prevalence as compared to gnomAD r. 2.1. RESULTS: An estimated 9.3% of discovered variants were not present in dbSNP. We report statistically significant overrepresentation of pathogenic variants for several Mendelian disorders, including phenylketonuria (PAH, rs5030858), Wilson's disease (ATP7B, rs76151636), factor VII deficiency (F7, rs36209567), kyphoscoliosis type of Ehlers-Danlos syndrome (FKBP14, rs542489955), and several other recessive pathologies. We also make primary estimates of monogenic disease incidence in the population, with retinal dystrophy, cystic fibrosis, and phenylketonuria being the most frequent AR pathologies. CONCLUSION: Our observations demonstrate the utility of population-specific allele frequency data to the diagnosis of monogenic disorders using high-throughput technologies.